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Abstract. We give explicit index-free formulae for all the degree six (and also 
degree four and two) algebraically independent local unitary invariant polynomials 
for finite dimensional fe-partite pure and mixed quantum states. We carry out 
this by the use of graph-technical methods, which provides illustrations for this 
abstract topic. 
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1. Introduction 

The notion of entanglement of a composite quantum system is known to be invariant 
under unitary transformations on the subsystems, so the investigation of local unitary 
(LU) invariants is a natural way of studying quantum entanglement. In this paper, 
we give illustrations for the general results of Hero et. al. [TJ [2] and Vrana [3) [4] on 
LU-invariant polynomials for pure quantum states. In [4], it has been pointed out that 
the inverse limit (in the local dimensions) of algebras of LU-invariant polynomials of 
finite dimensional fc-partite quantum systems is free, and an algebraically independent 
generating set for that has been given. This approach using the inverse limit 
construction is different from the usual, when the LU orbit structure is investigated 
first — for given local dimensions — and then invariants separating the orbits are being 
searched for The structure of algebras of LU-invariant polynomials 

for given local dimensions is very complicated, the inverse limit of these [U II], 
however, has a remarkably simple structure: it is free 0J, and an algebraically 
independent generating set can be given for that. Moreover, from the results for 
pure states, one can also obtain algebraically independent LU-invariant polynomials 
for mixed states [4]. 

The aim of this paper is to draw the attention of researchers working in the field of 
quantum information to the approach above — provided by researchers with expertise 
on representation theory — by hinting at the nature of results obtained within this 
approach. In particular, we write out explicitly the linearly independent basis of the 
inverse limit of algebras and single out the members of the algebraically independent 
generating set from them in the first three graded subspace of the algebras. We give 
these polynomials in an index-free form for arbitrary number of subsystems. 

The outline of this paper is as follows. In section we introduce the general 
writings of an LU-invariant polynomial and preclude the appearance of identical ones 
in a less abstract way than was done originally in [2 0]. We discuss the cases of 
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pure and mixed quantum states. In section [31 following [HH], we introduce graphs 
for the LU-invariant polynomials. Then we learn to read off matrix operations (such 
as partial trace, matrix product, tensorial product or partial transpose) from graphs. 
If this can be done for a whole graph of an LU-invariant polynomial, then we can 
write a nice index- free formula for that by these operations. In section [4] we give 
these index-free formulae for pure state invariants of degree two, four and six. Using 
graphs, these formulae can be given for arbitrary number of subsystems. In section[5j 
we show the formulae for mixed states and we discuss the connection of pure and 
mixed quantum states from another point of wiev. In section [51 we give an algorithm 
for the construction of the labelling of different invariant polynomials of degree six. 
(For degree two and four, this task is trivial.) Summary and some notes are left for 
section [Tj 

2. Local unitary invariant polynomials 

2.1. Invariants for pure states 

Let H = Hi <g) ... ® T-Lk be the Hilbert-space of a fc-partite composite system, where 
dim'Hj = nj and n denotes the fc-tuple of these local dimensions: n = (m, . . . ,«&)■ 
An element of the Hilbert-space can be written as \ip) = ipi lt ... t i k \ii, ■ ■ ■ ,ife), where 
\ij) G 7ij for ij — l,...,nj is an orthonormal basis for all 1 < j < fc, and the 
summation over ij = 1, ...,rij is understood. As usual in the topic of quantum 
invariants, the norm of i\) does not have to be fixed. 

It is well-known (see e.g. in [9]) that the way to get local unitary invariants is the 
following. We write down the term (V'ti,...,**'?^' z' ) m times (with different indices) 
and contract all primed indices with unprimed indices on the same Wj . A polynomial 
obtained in this way is of degree 2m, — degree m in the coefficients and also in their 
complex conjugates. This is the only case in which unitary invariants can arise [4], 
so it is convenient to use this natural gradation, and to call this polynomial of grade 
m. (In the case of mixed states the grade coincides with the degree in the matrix- 
elements of the density matrix.) The possible index-contractions on an Hj are encoded 
by the elements of S m , the group of the permutations of m letters, (j j G S m tells us 
that the primed index of the Zth term is contracted with the unprimed index of the 
(Tj(Z)'th term, so there is an index-contraction scheme for all fc-tuples of permutations 
a = (<7i, . . . ,<7 fc ) £ S^: 

fer{lp) =1pi],...,il •••V'irV.-,ir*^W 1) ,-- fc (D ' ' " ^ -i (™) ,"*(«') 5 (1) 

where the summation over — l,...,nj for all 1 < j < k and 1 < I < m is 
understood. (The lower labels of i's refer to the subsystems and the upper ones refer 
to the different index-contractions.) 

However, different fc-tuples of permutations can give rise to the same polynomial. 
We have the terms (i/y ip.axtf) ) 171 times: 

(A{,...i>%iM i .J(i } ii,...i> i °iw^) ■ ■ ■ fervi-ii-)^,): 

but it makes no difference if we permute the ipy 's or ^-"lO ' s among these terms, 
since, being scalar variables, they commute. This is equivalent to the relabelling of the 
indices (in the upper labels), which can be formulated by the permutations a, (3 G S m 
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encoding the permutations of tp.^i) 's and ipd 's, respectively: 

i 1 .... li--- 

(Al... V'-iw^. ,)(^,...^iffl v J ■ ■ • (^ir,...^-^...) 

= (^if 1) ,...^i c "' i(1) ,...)(^if (2) ,...^ij"'i w ,...) • • ■ (^if <m) ,..A°" i(m \...) 

= (fe,..>.» 1 fl- i (D )(^,...V , . QCT1 0-i ( 2) ) • • • (^....VL^-icm) )• 
*i >••• i i *i >••• 

(Here we have written out only the indices on Hi to get shorter expressions, but, 
obviously, the same a and /3 work on every index running on every Hj.) Therefore we 
have 

f{a u ...^ k ){^) = /(a<ri0- 1 ,...,e«Tk i 8- 1 )W')i ( 2 ) 

giving rise to an equivalence relation on S^: 

a ~ a' iff 3a, (3 S 5 m : cr'- = a<jj/3~ , 1 < j < (3) 

and the equivalence classes are denoted by [cr]^ — [ai, . . . , Uk]~ = 
{(acri/3 -1 , . . . , \ a, (3 £ S m }. The set of these equivalence classes is 

the double-cosets of by the diagonal action: A\5^/A, where the subgroup 
A = {(5,...,5)\5eS m }CS^. 

Thus, the ambiguity arising from the commutativity of the m terms 0... and 
"0... m ® has been handled by the labelling of the polynomials by the elements of 
A\S^/A. As a next step, it would be desirable to get one representing element 
for every equivalence class. Unfortunately, this can not be done generally, (i.e., for 
an arbitrary m,) but we can make the equivalence classes smaller by throwing off 
some of their elements in a general way. Every equivalence class has elements having 
the identity permutation e in the last position. Indeed, we have ao~uP~ l = e in 
(a<7i/3 -1 , . . . , if we set a — ficr^ 1 : 

(atri/3" 1 , . . . ,aa k -ifi~ x ,aa k f3~ x ) ~ 

{pa^<T X p-\ • • • , P<Tk e) = • • • , ^fc-i/?" 1 , e), 

which is actually an orbit of S*, - 1 x {e} under the action of simultaneous conjugation. 
So it is useful to define another equivalence relation on 5^ 

a « a' iff 3/3 g S m : c^ = /Soj-jS -1 , 1 < j < k (4) 

and the equivalence classes arc denoted by [a]~ = [cti, . . . , = 
{(/3cti/3 _1 , . . . , /3(Tfc/3 _1 ) | (3 e S m }. The set of these equivalence classes is denoted 
by Sm/S m . This equivalence is defined on S^ 1 in the same way. S*^ 1 can be in- 
jected into S!^ by i : S^ 1 c — > S 1 ^ as «(<Ti, . . . , Ofe-i) = (cti, . . . , (t^-i, e), which is 
compatible with the equivalence «, but not with ~. Since if a « cr' then cr ~ cr', a 
^-equivalence class is the union of disjoint ^-equivalence classes 

M~ = U [cr (2) ] K U . . . (5) 

The elements of [a] ~ which have Uk — e form one of the ^-equivalence classes of the 

right-hand side. This pa-equivalence class (element of S^ 1 / S m ) is also suitable for 

the labelling of the polynomials instead of the original ^-equivalence class (element 

ofA\5*/A). " 

The meaning of the choice ao~kf3 1 = e is that the indices on Hk are contracted 

inside every term (ipii j ip-°iW -i )■ This "couples together" the pairs of ip and ib. 

k l \ '■■■> l k 

The simultaneous conjugation means the permutation of the m terms {ip,..ip ), which 
is the remaining ambiguity arising from the commutativity of these terms. Note, that 
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we have singled out the last Hilbert-space Wk in this construction. In the general 
aspects, it makes no difference which Hilbert-space is singled out, but as we write 
the pure-state invariants using matrix operations, it can happen — and usually it will 
happen — that this freedom manifests itself in the different writings of the same pure 
state invariant. 

Summing up, for a composite system of k subsystems, the LU-invariant 
polynomial given by [a\, . . . , ak-i]& £ S 1 ^ 1 / S m is 

By the use of S^ 1 / S m labelling, we have got rid of the formal equivalence of 
polynomials arising from the commutativity of the terms, and have got a set of LU- 
invariant polynomial for the elements of the set S^ 1 /S m . Can it happen that different 
elements of S^ 1 / S m gives the same polynomial? Are there linear dependencies among 
these polynomials? It is not known in general, but sometimes there is more to be 
known: ((6]) gives a linearly independent basis in each to graded subspace of the inverse 
limit of the algebras [2J. Moreover, — as the main result of [3] states, — an algebraically 
independent generating set is formed by the polynomials given in (JS]) for which the 
defining k — 1 permutations together act transitively on the set of m labels. For the 
algebras of given local dimensions n — [n\ , . . . ,nk), the above polynomials form a basis 
as long as to < nj (for all j), otherwise they become linearly dependent. The algebraic 
independency also fails if we restrict ourselves to given local dimensions. (The algebra 
of LU-invariant polynomials is usually not even free for given local dimensions.) 



2.2. Invariants for mixed states 

Now consider a mixed quantum state of the fc-partite composite system. 
This state is given by the density operator, acting on Ti, written as g = 

Qii,...,i k ;i{,...,i'.\h, ■ • ■ ,ik)(i[, ■ • ■ ,i'k\- The density operator, by definition, a positive 
definite self adjoint operator, but, as usual in the topic of quantum invariants, the 
trace of g does not have to be fixed. 

The general form of an LU-invariant polynomial is given by a simillar index- 
contraction scheme, encoded by a = (ax, . . . , at) £ S^ n , as in the case of pure states: 

where the summation over ij = 1, . . . , rij for all 1 < j < k and 1 < I < to is 
understood. (We denote the pure and the mixed state invariants with the same symbol, 
the distinction between them is their arguments: they are vectors and matrices, 
respectively.) 

Here we can carry out a similar construction as in the case of pure states, with 
one difference: the building blocks of the polynomials are the (Qi 1) ... ! i k -i' ll ... t i l ) matrix- 
elements of the density operator instead of the former (ipi 1 ,...,i k ijj i ' i $/ )'s. Hence there 
is no step corresponding to the "double coset" construction: we can not move the "two 
parts" of g independently as has been done in the case of ijjtp, since in general g is not 
of rank one. This means that we can not relabel the primed and unprimed indices 
independently. The possible relabelling is given by /? £ S m : 

( ^ ' ^ ;....; ^ " l(1) ,...)( ^ '^2,...;^" 1<2, ,...) - ' ' ,...;i^ (m) ,..) 

= (^W^^^lW^JWifW^^^lW^J ■ ■ • (£j^ m ),... ; jf CT l< m ),...) 
— (£,-1 .^"if-Hl) )(Q,2 .-/3^i3- 1 (2) )...{Q p^p-Urn) )• 
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Therefore we have 

f(ai,...,<T k )(Q) = /( / 3<T 1 /3- 1 ,...,/3<T fe /3- 1 )(£') J (8) 

the elements of the orbits in under the action of simultaneous conjugation gives 
the same polynomial. Let these orbits be denoted by [o~i, . . . , ak]^, as before, and the 
LU-invariant polynomial given by this is 

f{<Tl = Q.l ,-1.-1(1) a k (l) ■ ■ ■ Q , m .,<ri(m) .<r k (m) . (9) 

The independency of these follows from the independency of the pure state 
invariants when that is the case for the latter ones. This is because we can 
obtain the independent mixed state invariants of the system with local dimensions 
n = (m, . . . ,nk), if we add a large enough Hk+i Hilbert-space, and calculate the 
invariants © for a pure state \<f>) £ H^Hk+i- (See [I] for the abstract construction.) 
Since in © we have not permuted the last (this time k + l'th) indices, we can read off 
the invarians for g = Tr^+i |^>)(<^| from ©. (If dimH^+i > J\j=i dirn'W,-, then g can 
be of full rank, and we can get all g acting on H in this way.) Note that if we simply 
substitute g by a pure state \tp)(ip\ in ©, then we do not get a linearly independent set 
of fc-partite pure state invariants for all the labels [<ti, . . . , <7fc] ra £ S!^/S m . However, 
if we restrict this for the case when o~k — e, then we get back the linearly independent 
set of pure state invariants from the linearly independent set of mixed state ones of a 
fc-partite system: 

/[«n,...,<*-i]« wo = f\.u-^-^]M)(M) = ^i,..,^_ I ]«c&* mm- (io) 

3. Graphs and matrix operations 

3.1. Graphs of invariants 

The index-contraction scheme of the LU-invariant polynomials given in the previous 

section can be made more expressive by the use of graphs [H [4]. For a grade m 

invariant, given by a — (ct 1: . . . , af.) S S 1 ^, one can draw a graph with m vertices with 

the labels 1 < I < m. These vertices represent the m terms (Vy ji ip-<rim -"kC)) m 

i>---i fc » x ,...,t fc 

|T]). The edges of the graph are directed and coloured with k different colours. The 
edges of the j'th colour encode the index contractions on the j'th Hilbert-space given 
by the permutation <jj of (a±, . . . , Uk)'- for every 1 < I < m there is an edge with head 
and tail on the Z'th and <7j(Z)'th vertex, respectively, meaning a contracted j'th index 
of the Z'th ip and the Gj(l)'th ip. 

As we have seen in the previous section, some elements of give rise to the 
same polynomials, as was elaborated in the previous section. How can we tell that 
story in the language of graphs? Following the previous section, for a given a € 5*,, 
the elements of the ^-equivalence class [a]~ G A\Sf n / A give the same invariant. First 
we set (Tfc = e, which means that we select the graphs having a loop of colour k on 
every vertex. For a given ^-equivalence class, there are still many graphs of that kind, 
and they are given by the elements of the corresponding ^-equivalence class. How 
are these graphs related to each other? The simultaneously conjugation by a f3 £ S m 
means the relabelling of the vertices, that is, the relabelling of the indices (in the 
upper label) of the terms (Vy ji i/ , "i(o "tC))- So the elements of a ^-equivalence 

class give the same graph with all the possible labellings, and the ^-equivalence class 
itself gives an unlabelled graph. 
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Since the elements of Sf n related by simultaneous conjugation gives rise to the 
same unlabelled graph, the decomposition in shows that there may exist many 
unlabelled graphs (many ^-classes) giving rise to the same polynomial defined by a 
given ^-class. For example, we can set o~j = e for a j ^ k, which results graphs where 
the edges of colour j ^ k form loops on every vertex. On the other hand, there may 
be ^-classes in the given ^-class which does not contain e. All of these graphs give 
the same polynomial, but it can happen that some of them can be formulated using 
matrix operations (in different ways for different graphs) and some of them not. (It 
turns out (see in next section) that every polynomial can be formulated using matrix 
operations up to m = 3.) 

The case of mixed states is simpler, because there are no ~-classes involved. 
For a grade m invariant given by a £ S^, the vertices represent the terms 
(g. t .^(j) .»».(!))) an d only the polynomials given by the elements of the w- 

equivalence class [cr]« € S^/ ' S m are the same by the commutativity of these terms. 
This means that we simply throw off the labelling of the vertices of the graph given 
by cr. 

3.2. Graphs of matrix operations 

The building blocks (V , u,...,i fc ' ! /'i' 1 ,...,^) an d {Qi 1 i k ;i' 1 ,...,i' h ) of the polynomials are 

matrices with row and column indices being the unprimed and primed ones, 
respectively. So we expect that some of the invariant polynomials can be written using 
only matrix operations, such as partial trace, matrix products, tensorial products 
or partial transpose. How can we read off matrix operations from the graphs 
corresponding to the invariant polynomials? This is a difficult question in general 
and — as we will see — not all graphs can be encoded using matrix operations. It is 
more instructive to look at the graphs coming from the matrix operations first, and 
then to search for these elementary subgraphs in a general graph coming from a 
polynomial given by an element of 

Let us see some matrix operations and their graphs. The matrix multiplication 
means contraction of the column indices of the first matrix with the row indices of 
the second matrix, the trace means contraction of the column indices with the row 
indices and the partial transposition means the swap of the given row and column 
indices. First consider only the indices belonging to the Hilbert space of only the first 
subsystem — i.e., we have edges of only one colour. For a general matrix M, which is 
represented by the vertices of the graphs, the multiplicaton by itself gives the edge 
from one vertex to another, the r'th power M r is a chain of edges (without loops), 
and the trace of it closes this chain into a loop. (See in the first row of figure [TJ) 
Now let us take into account indices belonging to the second subsystem. TrM r 
is the same loop as before but with doubled edges, while the partial transposition 
Tr(A/ T2 ) r reverses the loop of the corresponding colour (second row of figure [T]). 
The partial traces in TrTri(M r ) Tri(M s ) make smaller loops on a subsystem (third 
row of figure [l}. There is a little trick, which is prooved to be very useful later: 
Tr(Tri M 2 ) (Tr : M) = Tr M 2 (It ® Tri M). On the language of graphs we just bend the 
corresponding edge next to the vertex representing Tri M, and we draw a circle on it, 
representing the identity matrix, which is just contract indices (last row of figure [T]). 
If the graph is the union of disjoint graphs, then the corresponding polynomial is 
factorizable, since the summations corresponding to the disjoint pieces can be carried 
out independently. This almost trivial situation is getting more complicated, if we take 
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into account the indices of all the subsystems — i.e., the edges of all colours. Examples 
are shown in the next section. 



M 2 : 



TrM 2 : 



TrM r : 
1 2 




Tr(M T2 ) r : 
1 2 3 



Tr(Tri M r )(Tri M s ): 
12 r r 




Tr(Tri M 2 )(Tri M) 



TrM 2 (Ii <g> Tri M) 



'•O 



Figure 1. Elementary matrix operations represented by graphs. In the first 
row: k = 1, there is only one colour of edges representing the index contractions. 
In the other rows: k = 2, two different colours of edges correspond to the index 
contractions on the two Hilbert-spaces: black and red on Hi and H2, respectively. 
(Black and dark grey in greyscale printing.) 



4. Pure state invariants 

In the following, we illustrate how a pure state LU-invariant polynomial (encoded by 
[a]~ £ A\S* /A) is given by different unlabelled graphs (encoded by [a}~ € 5^/5 m ). 
(While an unlabelled graph is given by different labelled graphs (index-contraction 
scheme, encoded by a £ 5^)0 O n the other hand, different unlabelled graphs give rise 
to different writings by matrix operations of the same polynomial. The polynomials are 
labelled here by the elements of S*^ 1 /S m instead of the elements of A\5^/A, which 
give special unlabelled graphs having loops of colour k on every vertex. (Sometimes, 
e.g., in [2], these loops are omitted, and only the first k — 1-coloured edges are 



drawn.) For a permutation k — 1-tuple a £ 
write for the corresponding invariant [i(cr) 



ok- 

U = [t(<r)U U [a^} x U [cr( 3 )]„ U 



[a]~ £ S 1 ^ 1 /S m , and we can 



where 



£ iS^ are representing elements of ^-classes giving different graphs for 
the same invariant. 

Let us see how these technics work. As a warm-up, we show for all k the trivial 
case of in = 1 and the almost trivial case of m = 2. This is followed by the case 



All degree six local unitary invariants of k qudits 



8 



of m = 3, which is more interesting because of the non-Abelian structure of S3. 
This is done for all k too. Let \ip) G "H, and, as we have seen, everything can be 
formulated using the rank-one projector tt = 7r 12 ...fe = We denote the reduced 

density matrices with the label of subsystems which are not traced out, for example 
7T2...fc = Tri 7Ti 2 ...fc etc. 

Invariant polynomials of grade m = 1 (degree two) 

For m = 1, we have the trivial S\ = {e}, and for all k number of subsystems 
[e, e, ...,e]^ = [e,e, ...,e]~ (A\Sf/A = S^/Si = Si) meaning only one kind 
of graphs, having only one vertex. Every edge — of k different colors for the k 
subsystems — starts and ends here. (See in figure [2]) This graphs mean a simple 
trace, which is the only possible index contraction. The label of the polynomial is the 
only one element of S^ 1 /Si = Si, and 

/[e,..., e]si W=Tr^ 2 ... fe = |^|| 2 . (11) 

k = 1: k = 2: k = 3: arbitrary k: 

[e]w [e,e]~ [e,e,e]~ [e, e, e, . . . , e]~ 

Q# % 

Figure 2. Graphs corresponding to the m = 1 invariant polynomials. Black, 
red, blue and green edges (from the darkest to the lightest in greyscale printing) 
represent index-contractions on the first, second, third and last Hilbert spaces, 
respectively. 



4-2. Invariant polynomials of grade m = 2 (degree four) 

For m — 2, we have £2 = {e,i} (where e = (1)(2) and t = (12)) with the conjugacy- 
classes [e] and [t], so the labels of the polynomials are S^ -1 / S2 — S^ 1 for all k. On 
the other hand, [a]^ = [a)~ U [ctJw, (where <fj = 57, and ? = e, e = t) so there are two 
kinds of graphs for every polynomial. 

For one-partite system, (k = 1, iz = 7ri,) the only polynomial is for 

[e]~ = [e]*U[t]«. 

From its graphs, (see in figure [31) we have 

/ []si (V) = (Tr7r 1 ) 2 =Tr7r 1 = W 4 . 

For two-partite system, (fc = 2, 7r = 7Ti2,) there are two linearly independent 
polynomials. These are given by 

[e ) e]^ = [e,e]«U[t,t]« ) 
[i,e]^ = [t,e]« U [e,t]«. 

From their graphs, (see in figure [3j) we have 

/ [e]K W = (Tr7r 12 ) 2 = Tr7r? 2 ^||^|| 4 , 
/[t]«W0 = TrTT 2 = TrTr 2 . 
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For three-partite system, (fc = 3, 7r = 7r 12 3,) there are four linearly independent 
polynomials. These are given by 

[e,e,e]^ = [e,e,e]~ U [t,t,t}~, 

[e, t, e]~ = [e, t, e]~ U [t, e, t]~, 

[t, e, e]^ = [t, e, e]~ U [e, i, t}~, 

[t,t,e]„ = [t,*,e]„ U [e.e.t^. 

From their graphs, we have 

/[e,e]„ WO - (Tr7T 123 ) 2 = It 7T? a3 = ||V|| 4 , 

f[e,t] a {4>) = Tr 7r 2 = Tr7r 2 3 , 

/[t.e^O) = Tr7r l = Tr7r 23> 

/[t,*] M W = Tr7r 2 2 = TrTr 2 . 



fc = 1: 







(Tr^) 2 = TrTr 2 

fc = 2: 

[e,e]„ ~ [t,i]„ 

• • • ^ ^ • 

(TrTr) 2 = TrTr 2 

M« ~ [e,t]« 

Tr(Tr 2 tt) 2 = Tr(Tri tt) 2 

Figure 3. Graphs corresponding to the m = 2 invariant polynomials for k = 1 
and 2. Black and red edges (black and dark grey in greyscalc printing) represent 
index-contractions on the first and second Hilbert spaces, respectively. 

The construction of these formulae can be easily generalized to arbitrary number 
of subsystems. For this, take a look at the left graph of the last line of figure [31 This 
time, let the red lines (dark grey in black and white printing) represent the index- 
contractions on all Hilbert-spaces on which <Tj = e, and the black lines represent the 
index-contractions on all Hilbert-spaces on which <jj = t. Thus, we can read off the 
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matrix operations for arbitrary k. The other way of writing the polynomial can be 
reached by the interchange of the roles of the black and red lines. So, for arbitrary 
number of subsystems (tt = 7r 12 ...fe), for the polynomials for [cr]^ = [cr]~ U [a}~, we 
have 

f[a 1 ,....a k ^ 1 ] a W = Tr(Tr {fc}u{iK=e} tt) 2 = Tr(Tr {iki=t} tt) 2 . (12) 

This was well-known for qubits [9] . The number of these — the dimension of the grade 
rn = 2 subspace of the inverse limit of the algebras — is 2 k ~ 1 . The set of algebraically 
independent generators contains all the m = 2 polynomials from (|12j) . except the ones 
for which there are only e's in [a%, . . . , (7fc_i]~ labelling the polynomial. (This is the 
only way for the permutations not to act transitively on the set of m — 2 labels.) The 
number of these is 2 fc ~ 1 — 1. 

4-3. Invariant polynomials of grade m = 3 (degree six) 

For 77i = 3, we have S 3 = {e, s, s 2 1 t,ts,ts 2 } (where e = (1)(2)(3), t = (12)(3), and 
s = (123)) with the conjugacy-classes [e] = {e}, [s] — {s,s 2 } and [t] — {t,ts,ts 2 }. 
This time, we have no simple general rule for the splitting of a ~-class to ^-classes. 
For one-partite system, (k = 1, tt = 7Ti,) the only polynomial is for 

[e]« = [e]„U[t]«U[i]«. 

From its graphs, (see in figure 01) we have 

/D»(V0 = (Tr 7Ti) 3 = Tr tt 2 Tr tt\ — Tr tt\ = ||t/-|| 6 . 

For two-partite system, (k — 2, tt = tt\2,) there are three linearly independent 
polynomials. These are given by 

[e,e]„ = [e,e]«U [t,t]« U [«,*]«, 

[t, e]„ = [t, e]~ U [e, t}~ U [s, t]~ U [t, s]~, 

[s, e]^ = [s, e]~ U [e, s]» U [s, s 2 }~ U [t, is]~ 

(so S3/S3 ~ {[e], [s], [i]}). From their graphs, (see in figure IH) we have 

f [e] JiP) = (Tr^ 12 ) 3 = Tr7r 2 2 Tr^ 12 = Tmf 2 = |^f, 

/M»W0 = (Tr7r a ) 2 Tr7Ti2 = Tr(Tr a 7r 2 2 )7T6, 

/Mr.WO = Tr7r a = TrC^S ) 3 = ^(^i ® ^2)^12, 
for all a, b 6 {1, 2}, a ^ 6. 

For three-partite system, (k = 3, 7r = 7i"i23,) it turns out that there are eleven 
linearly independent polynomials. These are given by 

[e, e, e]^ = [e, e, e]~ U [t, t, t]~ U [s, s, s]~, 

[e, i, e]^ = [e, rj, e]~ U [i, e, t]~ U [s, £, s]~ U [t, s, /;]~, 

[t, e, e]^ = [<, e, e]~ U [e, <, t]~ U [t, s, s]~ U [s, t, 7;]~, 

[i,t,e]~ = [t,t,e}~ U [e,e,t](a U [s,s,i]~ U [*,*,«]«, 

[e, s, e]^ = [e, s, e]~ U [s, e, s]~ U [s, s 2 , s]~ U [i, ts, t]~, 

[s, e, e]^ = [s, e, e]~ U [e, s, s]~ U [s 2 , s, s]~ U [ts, t, t]~, 

[s, s, = [s, s, e]~ U [e, e, s]~ U [s, s, s 2 ]~ U [t, i, tsj^, 

[s, s 2 , e]^ = [s, s 2 , e]~ U [s, e, s 2 ]~ U [e, s, s 2 ]~ U [t, is, £s 2 ]~, 

[i, s, e]^ = [i, s, e]~ U [t, e, s]~ U [e, t, is]~ U [t, s, s 2 ]~ U [s, t, ts}~ U [s, t, 
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[e]„ ~ [t]„ ~ [sU 

• • • 

(TVtt) 3 = TW 2 Trvr = TVtt 3 
ft = 2: 

[e,e]„ ~ [f,t]« ~ [s,s}~ 

5 A /A 

• • • ^ ^ • • • 

(Tr^) 3 = TW 2 Tr7r = TVtt 3 

M» ~ [e,t]„ - [*,*]» ~ [t,a]« 



<v) <v) 

Tr(Tr 2 7r) 2 Trvr = Tr(Tri tt) 2 Tr vr = Tr(Tr 2 tt 2 Tr 2 tt) = Tr(Tri tt 2 Tri vr) 

[s, e]~ ~ [e, s]~ ~ [s,s 2 ]ps ~ [i, is]fs 





^ — " ^5 ^ — ' ^ ^ — 

Tr(Tr 2 7r) 3 = Tr(T ri tt) 3 = Tr(7r T2 ) 3 = Tr(I a <g> Tri 7r)7r(Tr 2 tt ® I 2 ) 

Figure 4. Graphs corresponding to the m = 3 invariant polynomials for 
k = 1 and 2. Black and red edges (black and dark grey in greyscale printing) 
represent index-contractions on the first and second Hilbert spaces, respectively. 
The formulae of the polynomials given by matrix operations, which can be red off 
from the graphs, are also written out. For the last one, we have use the trick in 
the last line of figure [l] twice. 
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[s, t, e]„ = [s, t, e}~ U [e, t, s}~ U [t, e, ts}~ U [s, t, s 2 ]~ U [t, s, ts]~ U [i, s, ts 
[t, ts, e]^ — [t, ts, e}~ U [e, s, i]~ U [s, e, i]~ U [s, s 2 , t]~ U [t, ts, s]~ U [t, ts 2 , 

Here we do not write out all the 49 formulae for the graphs coming from the ^-classes 
above, we just show two or three of them for every polynomial. 



/[e,e]„0) 


= H^ll 6 = 


: (Tr7Ti2 3 ) 3 , 








= Tr 7T123 


Tr 7r 2 = Tr7ri2 3 Tr7r 2 3 , 








= Tr 7T123 


Tr 7r 2 ee Tr 7Ti 2 3 Tr 7723, 






f[t,t]jip) 


= Tr 7r 123 


Tr7r 12 = Tr7r 123Tl'7r 2 , 






/[e,s]«W>) 


= Tr7r 2 3 E 








/[3,e] K (V') 


= Tr = 


= Trvr 2 3 3, 






/[«,«]« WO 


= Trwl 2 : 


EETrvr 3 , 






/[*,*]« W) 


= Tr(h « 


)7r 2 )7r 2 2 ee Tr(Ii ® 7r 3 )7r 2 3 ee 


Tr(7r 2 $ 


5 ^3)^23) 


W) 


= Tr(7n $ 


5 I 2 )7T 2 2 EE Tr(I 2 ® 7r 3 )7T^ 3 EE 


Tr(7ri (> 


5 7T3)7Tl3, 


/[t,t»]« WO 


= Tr(7n $ 


5 7r 2 )7ri2 ee Tr(7Ti ® I 3 )7r 2 3 ee 


■ Tr(7T 2 C 


g>l3)7T23, 


/[«,a a ]«W0 


= Tr(^ 


) 3 a, b e {1,2,3} distinct. 







The last one of these is the Kempe- invariant [T2], which has arisen in the context 
of hidden nonlocality. Kempe has defined this for n = (2,2,2), i.e., three qubits, 
but it can be written by her definition (see (17) in [12]) for all n — (711,713,713) 
three qudit systems. It is observed by Sudbery [5] that for three qubits the Kempe- 
invariant can be expressed as f[ s , s 2 ]sl (ip) = 3/[ M]ai (V0 - /[ e , s ]„ (VO - f[s, s ] a (^) = 
3 /[«,t]«W0 - f[s,e] K {^) ~ f[s,suW = 3/[t >ts ] K (V0 - f[e,s] a (ip) ~ /[a, e ]«(^)- However, 
this is only for qubits: if m < for all j (so at least qutrits) then the 11 
polynomials listed above are linearly independent. Another important three-qubit 
permutation- and LU-invariant polynomial of degree six, which has arisen in twistor- 
geometric [T3] and Freudenthal [TJ] approach of three-qubit entanglement, is the 
norm square of the Freudenthal-dual of tp. It can be written as 6\\T(tp,ip,ip)\\ 2 = 
4f[ s ^u{ip) + 5/[ e , e ] K (-0) - 3/[ M]s3 (V>) - 3/[t ie]a! (V>) - 3f[t,t] a (tp)- Note, that this 
expression is not unique, since these /u^r polynomials are not linearly independent 
in the case of three qubits. 

It is not obvious, but the construction of the grade m = 3 polynomials can be 
generalized to arbitrary number of subsystems. To do this, consider an invariant 
given by o~ £ S3 , where o~j G {t, ts, ts 2 }. This can be seen in figure[51 with the evident 
redefinitions of the meaning of the colours: let the black, red, and blue edges (black, 
dark grey and light grey in greyscale printing) represent the index-contractions on all 
Hilbert-spaces on which crj — t, ts and ts 2 , respectively. Using the trick in the last 
line of figure [T] three times, we have, that the polynomial is given by the trace of the 
product of the three factors: I{j|a-,-=t} ® T^{j\ a:j =t} T) and other two with ts and ts 2 . 
Note, that the order of these are arbitrary, since it is related to the relabelling of the 
vertices of the graph. However, if there are some subsystems on which Uj = s, that 
fixes the order of these terms up to cyclic permutation. It turns out that we have to 
use the reverse ordering in the product, the terms liju.—^-y (g)Tr/ J -| .._ T i ir with r = ts 2 
first, then with r = ts and then with r = t, since they have the fixed point 1, 2 and 
3, respectively. On the subsystems on which o-j — s, the indices are intact, — partial 
traces act only on other subsystems, — they are contracted in the appropriate way. If 
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t,ts,ts 2 




t, ts, ts 2 








Figure 5. Graphs corresponding to the m = 3 invariant polynomials. Black, 
red, blue and green edges (from the darkest to the lightest in greyscale printing) 
represent index-contractions on the Hilbert spaces on which aj = t, ts, ts 2 and s 
respectively. For the first graph, we show how the trick in the last line of figure [l] 
was used three times. 



there are some subsystems on which <7j = s 2 , then we use tt { J i CT 3= s2 > instead of tt 
to reduce the situation to the known case. Similarly, if there are some subsystems 
on which o~j — e, then we use r Trr J -| .._ e i tt instead of tt. Summing up, for arbitrary 
number of subsystems (tt = ir\2...k) we have the following formula for the m = 3 
polynomials: 



T=tS 2 ,tS,t 

where the Y[ product symbol means non-commutative product, in the order of its 
subscript. This gives back the formulae for the special cases k = 1, 2 and 3. 

5. Mixed-state invariants 

In section [21 we considered the mixed state invariants of k subsystems as pure state 
invariants of k + 1 subsystems. By considering the graphs of the invariants, given in 
section [31 we can clarify this from another point of wiev. 

An invariant is given by a € S^ n : this encodes an index-contraction for the 
matrices of the operators tt = or g for pure or mixed states, respectively. If 

a 96 a' for a, a' € S 1 ^, then they give rise to different polynomials for mixed states, 
while it can happen, that a ~ a 1 , so they give rise to the same polynomial for pure 
states. In this case, the unlabelled graphs given by [cr]« and [o"']« ^ are related 
to each other by the independent permutation of the heads and tails of the edges, 
while the corresponding operation is the independent permutation of the coefficients 
ip and ifi in ([T]). This operation is not allowed for mixed states. In section [4] we have 
given the decompositions of ^-equivalence classes into ^-equivalence classes for some 
grade m and for some k numbers of subsystems, leading to the different writings of 
the same polynomial. For mixed states of k subsystems, these polynomials are not 
the same anymore. This offers us a different point of wiev, which seems to be more 
natural: let us consider the pure state invariants as the special cases of the mixed state 
invariants instead of considering the mixed state invariants as pure state invariants of 
a bigger system. We have the mixed state formula © for the set of invariants, if we 
substitute a pure state IV'XV'I hito them, then some of them will coincide, but we can 
keep this in hand. 

For the sake of completeness, we show the mixed state polynomials (j9)) below. 
Comparing these formulae with the ones for pure states, one can see how the fc-partite 
mixed state invariants are related to the k + 1-partite pure state ones (fTOj) . or, how the 




(13) 
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different fc-partite mixed state invariants coincide with the same fc-partite pure state 
invariants. 

Let g = Q12 t be a density matrix on H. The polynomials are labelled by 

5.1. Invariant polynomials of grade m = 1 (degree one) 
For m = 1, we have for all k 

fle,....e] a (Q) = Tr£>12...fc. (14) 

5.2. Invariant polynomials of grade m = 2 (degree two) 
For m = 2, for one-partite system, (fe = 1, g = g±,) we have 

/[e]„(e) = (Trgi) 2 , 

(for a ni = 2 one-qubit system, the determinant is an element of this subspace: 
2 det g = /| e L (p) — (p)) for two-partite system, (k = 2, g = g\2,) we have 

/[e,e]„(e) = (Trgi 2 ) 2 , 

/[e,i]»(e) = Tr^i 
/[t,e]„(e) = Trgi, 
/[*,*]„ (e) =Tr(? 2 2 , 
for three-partite system, (fc = 3, g = 0123,) we have 

f[e,e.e] a (g) = (Trgi 23 ) 2 , 

/[e.e.tJs.Ce) = Tr^i 

/[e,t,e]»(0) = Tr e2> 

/[e,i,i]«(e) = TreL 

/[t,e,e]« (e) = 1*01, 
/[t,e,t]„(e) = Tr 0i3, 
f[t,t,e]„(e) = Tr 012> 

/[t,t,t]«(e) = Tre 2 23 , 
and for arbitrary number of subsystems, (g = g\2...k) we have 

/[<ri,...,<7*]«(e) = Tr(Tr {ikj=e} g) 2 . (15) 
The number of these — the dimension of the grade m — 2 subspace of the inverse 
limit of algebras — is 2 fc . The set of algebraically independent generators contains all 
the m — 2 polynomials from (|15[) . except the ones for which there are only e's in 
[cti, . . . , <Tfc]ps labelling the polynomial. The number of these is 2 k — 1. 

5.3. Invariant polynomials of grade m = 3 (degree three) 

For m = 3, for one-partite system, (k = 1, £> = gi,) we have 
/ [e]=s ( e ) = (Tr ei ) 3 , 
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(for a ri\ = 3 one-qutrit system, the determinant is an element of this subspace: 
6detg = f[ e ] a (g) - 3f[ t ] x (g) + 2f[ s]x (g)) for two-partite system, (k = 2, g = g 12l ) we 
have 



fr i (n) 
/[e,e]« \y) 


— fTr nn"l 3 


J[e,t] a \Q) 


= Tr £>i2 Tr g\ , 


J[t,e]~(e) 


= Tr gi2 Trg x , 


/[t,t]«(e) 


= Trg 12 Tr^ 2 , 


/[e,«]«(e) 


= Tr4 


/[«,e]„(e) 


= Trg 3 , 


/[«,«]« (e) 


= Tr e ? 2 , 


/[*,«]« (e) 


= Tr(Ii ® g 2 )g 2 12 , 


/[«,*]» (e) 




/[t,t»]»(e) 


= Tr(£>i (g> £2)012, 


/[s,s 2 ] K (^) 


= Tr(^) 3 ^Tr(^) 3 



and for arbitrary number of subsystems (p = gi 2 ...k) we have 

/[ CT1) ... )CTfe ]«(e)= Tr J] (%k.=r}8)Tr {jV . e{T>e}} /«i^- 2 >),(16) 

where the J} product symbol means non-commutative product, in the order of its 
subscript. This gives back the formulae for the special cases k = 1 and 2. 

6. Algorithm for SyS 3 

The formula in (|13p gives grade m = 3 invariant polynomials for a (<ri, . . . , o~k-i) € 
5g _1 fc — 1-tuple of permutations, but the linearly independent ones are labelled by 
[<7i, . . . , <Tfe_i]~ e 1S3 l /Sa. Since the group-structure of S3 is not too complicated, we 
can give an algorithm to construct exactly one representative element a for all orbits 
[<r]~, i.e., to construct the elements of S3/S3. The choice r — k — 1 and r = k gives 
the labels for pure and mixed state invariants, respectively. 

Again, S3 = {e, s, s 2 , t, ts, ts 2 }, s = (123), t = (12)(3), and its conjugacy classes 
are [e] — {e}, [s] = {s, s 2 }, [t] — {t,ts,ts 2 }. First, we write the conjugation table: 



PIP' 1 ■ 



A 7 


e 


s 


s 2 


t 


ts 


ts 2 


e 


e 


s 


s 2 


t 


ts 


ts' 


s 


e 


s 


s 2 


ts 


ts 2 


t 


s 2 


e 


s 


s 2 


ts 2 


t 


ts 


t 


e 


s 2 


s 


t 


ts 2 


ts 


ts 


e 


s 2 


s 


ts 2 


ts 


t 


ts 2 


e 


s 2 


s 


ts 


t 


ts 2 



For every position of the list (oi, . . . , ay) g S3 there is a conjugacy class [<jj] of S3, 
which remains unchanged under simultaneous conjugation. For a given (a±, . . . , ay), 
we would like to single out one representative element in the orbit of simultaneous 
conjugation. To do this, we examine the orbits. 
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• If Oj = e for all j, we have the trivial orbit of length 1. Besides this case, we 
do not have to deal with the positions in which e's occur, since aj = e remains 
unchanged under simultaneous conjugation. 

• If (Tj G [s] (besides e) for all j, then we can choose the element of [a]~ which has 
s in the first position in which an element of [s] occurs. These orbits are of length 
2. 

• If Oj £ [t] (besides e) for all j, then we have two kinds of orbits. If aj is the same 
for all j for which <jj 6 [t], then we can choose the element which has t in the 
first position in which an element of [t] occurs. These orbits are of length 3. On 
the other hand, if there are at least two different Uj for which cr, £ [t] , then it is 
not enough to fix only one position. It can be checked by the conjugation table 
above that for the ordered pairs of different elements of [t] there exists exactly 
one permutation which brings them into (t, ts) by simultaneous conjugation. So 
we can uniquely choose the elements which have t in the first position in which 
an element of [t] occurs and ts in the first position in which a different element 
of [t] occurs. These orbits are of length 6. 

• If both [s] and [t] occur (besides e), then we have to fix two positions again. It 
can be checked by the conjugation table above that for every pair given by the 
elements of [s] x [t] there exists exactly one permutation which brings it into (s, t) 
by simultaneous conjugation. So we can uniquely choose the element which has 
s in the first position in which an element of [s] occurs and t in the first position 
in which an element of [t] occurs. These orbits are of length 6. 

With the help of the observations above, we can formulate the following algorithm 
generating S3/S3, i.e., the labels of the polynomials. 

(i) For every position of the list (01, . . . , oy) € 5J, assign one of the conjugacy-classes 
of 53. Do this in all possible ways, and apply the following steps for all of them. 

(ii) Write e into all positions to which [e] has been assigned. 

(iii) Take the first of the positions to which [s] has been assigned, and write s there. 
To the others of such positions, write either s or s 2 in all possible ways. 

(iv) If there is no position with [s], then take the first of the positions to which [t] 
has been assigned, and write t there. To the following of such positions, write 
cither t or ts in all possible ways, but after the occurrence of the first ts, write 
cither t, ts or ts 2 in all possible ways. On the other hand, if there is at least 
one position with [s], then take the first of the positions to which [t] has been 
assigned, and write t there. To the others of such positions, write either t, ts or 
ts in all possible ways. 

What is the number of the labels obtained in this way? This could be find by 
the use of some combinatorics, but we do not have to follow that way. If the local 
dimensions 3 < rij, then the elements of S^ -1 / S3 label the linearly independent grade 
m = 3 invariants, and their number, the dimension of the grade m = 3 subspace of 
the inverse limit of the algebras is given in [TJH]. For m = 3 pure state invariants, this 
is |S| -1 /5 3 | = 6 fe - 2 + 3 fc - 2 + 2 fe - 2 , for mixed states |Sf/S 3 | = 6 fe_1 +3 fe_1 + 2 fc_1 (see 
also in [15]). One can easily check that the set of algebraically independent generators 
contains all the m = 3 polynomials from ()13[) or (|16|) . except the ones for which — 
using the labelling algorithm above — there are only e's and t's in [<r]~ labelling the 
polynomial. (This is the way for the permutations not to act transitively on the set of 
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m = 3 labels.) The number of these is 6 k ~ 2 + 3 fc ~ 2 + 2 k - 2 - 2 k ~ 1 = 6 fc " 2 + 3 fc ~ 2 - 2 k -' 2 
for pure states, and 6 fc_1 + 3 fc_1 — 2 fe_1 for mixed states. 

7. Summary and notes 

In this paper we have written out explicitly the LU-invariant polynomials for pure and 
mixed states, given in ([6]) and ([9]), for to = 1, 2, 3. This was done for arbitrary number 
of subsystems of arbitrary dimensions. The key point — and the new feature — here is 
the independency [21 H] : the polynomials in ([5]) and © give the linearly independent 
basis of the to graded subspace of the algebra of LU-invariant polynomials if to < rij 
(for all j) [2] , and some of them — the ones for which the defining permutations together 
act transitively — become an algebraically independent generating set in the inverse 
limit of algebras, i.e., if nj — > oo for all j pQ. This independency result shows the 
power of the elegant approach using the inverse limit of the algebras of LU-invariant 
polynomials. 

However, for a given n = (m, ...,nu) system, it seems to be usual [6j [7] that 
it is not enough to use only the polynomials of maximal degree 2m, where m < rij 
for all j, for the separation of the LU-orbits. (According to the relatively simple case 
of n — (2,2,2) three qubits, where it is known [9) [TT], that we need an m = 3, an 
m = 4, and an m — 6 invariant polynomial — the Kempe invariant, the three-tangle, 
and the Grassl-invariant respectively, — beyond the m < 2 ones.) If to ^ rij for a j, 
the generators given in (|6]) and (|9]) will not be linearly independent, and the algebraic 
relations between them exhibit a complicated structure. 

New results are given by the nice compact formulae of grade m — 3 invariants 
for pure and mixed quantum states, (equations (|13p and (JTHJ), respectively) and in 
the algorithm generating the different equivalence classes of permutation fc-tuples of 
S*3 under simultaneous conjugation, given in section [6l The latter is necessary to 
eliminate identical polynomials. Connections between pure and mixed state invariant 
polynomials has been illustrated as well. These results are obtained by the use of 
graphs corresponding to the polynomials [21 0] . 

Note, that (i) the same degree of the pure state invariants in the coefficients 
and in their complex conjugate, (ii) the much simpler labelling of the mixed state 
invariants than the pure ones, (iii) considering the pure state invariants as the special 
cases of the mixed ones, seem to stress that the mixed states (density matrices) are 
the natural objects in the topic of unitary invariants instead of the pure states (state 
vectors). This approach is widely supported by the whole quantum physics, where the 
elements of the lattice of the subspaces of the Hilbert space are regarded to be more 
fundamental than the elements of the Hilbert space themselves. 

The illustrating polynomials given in this paper could have been written in a 
convenient form using partial trace, matrix product, tensorial product and partial 
transpose for grade m = 1, 2 and 3. However, we note that it can happen that a grade 
to > 4 invariant polynomial can not be written by using these operations only. At 
this time, we can not formulate general necessary and sufficient conditions for this, 
but we can give an enlightening example. For the use of matrix operations, we have 
to write down the matrices one after the other, this fixes the order of the vertices in 
some sense. The partial traces form loops of edges. If we can find an ordering of the 
vertices (up to cyclical permutations), in which these loops of every colours contain 
only adjacent points (with respect to this ordering) then the matrix-operations can 
be written for the entire polynomial. This situation can be seen in the third row of 



All degree six local unitary invariants of k audits 



18 
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c 




Figure 6. An example for an m = 4, k = 2 mixed state LU-invariant polynomial 
which can not be written by the considered simple matrix operations. 



figure Q] After some drawing, one can check that there is no such an ordering of the 
vertices for the graph in figure [51 which seems to be the most simple exapmle for such 
a situation. 
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